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We considered a spin chain with nearest neighbor and next nearest neighbor exchange interactions, 
anisotropic exchange interaction and Dzyaloshinskii-Moriya interaction. The conditions of the spiral 
spin state as the ground state were analyzed. Our purpose is to build the connection between the 
spiral state and the fully polarized state with a unitary transformation. Under this transformation, 
anisotropic exchange interaction and Dzyaloshinskii-Moriya interaction can be transformed to each 
other. Then we use positive semi-definite matrix theorem to identify the region of fully polarized 
state as the ground state for the transformed Hamiltonian, and it is the region of spiral spin state as 
the ground state of the original Hamiltonian. We also found that the effect of Dzyaloshinskii-Moriya 
interaction is important. Its strength is related to the pitch angle of spiral spins. Our method can 
be applied to coupled spin chains and two dimensional triangular lattice systems. Our results can 
be compared with the experiment data. 



I. INTRODUCTION 

The spiral spin states has been the subject of study 
for more than sixty years 1 - - — . Yet it still gives rise to sur- 
prising physical properties. We focus on the multiferroic 
phenomenon found recently in numerous compounds 4 . 
Experiments^ showed that in the multiferroic material 
the magnetic and ferroelectric orders are closely related. 
What is more intriguing is that only certain types of mag- 
netic orders, can couple to ferroelectricity. Furthermore, 
only spiral spins configuration gives rise to strong cou- 
pling between electric polarization and magnetic order. 
The spin-current model£ provided a plausible explana- 
tion for this phenomenon. 

Beside multiferroics, the spiral spin state was found 
in many other transition metal compounds. In the 
early study®' 9 , they have presented the evidences of 
the existence of spiral spin state in multiferroics. In 
LiC^O^ - — and NaC^Oz^, there are one-dimensional 
spin chains consist of edge-sharing octahedra. The Cu-O- 
Cu bond angle is almost 90-degree. This renders the su- 
perexchange interaction between nearest neighbor (NN) 
weak and ferromagnetic^ and the exchange between next 
nearest neighbor (NNN) not negligible. The spiral spin 
configuration can also be found in higher dimension sys- 
tems. The structure of ACr0 2 (A=Cu, Ag, Li, or Na)^ 
is a two dimensional triangular lattice, and its bond an- 
gle of Cu-O-Cu is also close to 90-degree. This can also 
be the cause of the spiral spin state. 

The existence of spiral spin configuration can be at- 
tributed to the frustration in the system. A relatively 
simple case is a one-dimensional spin chain with NN and 
NNN exchange interactions, to be denoted as J\ and J2 
respectively. Frustration is caused by their competing 
tendencies of aligning spins. This kind of systems is often 
called zigzag spin chains. There have been much analy- 
sis on this subject. Exact solutions have been found for 
special cases. The most notable case is the dimer state 
at J1/J2 = 2 found by Majumdar and Ghosh™,. At the 
other end J1/J2 = —4, it was found that the fully polar- 
ized (FP) state and uniformly distributed resonant va- 



lence bond stated are degenerate ground states (GS). In 
general, the phase diagram is summarized by Bursill 18 . 
The boundary of the frustrated region is identified by 
numerical calculation. White and Affleck— calculated 
the correlation function and provided solid ev- 

idence of the existence of the spiral spin state. Further- 
more, it has been found that in zigzag spin chains, there 
exists chiral order— ~—. This is another indication of the 
extensive existence of spiral spin state. Hence the exis- 
tence of spiral spin state becomes an important subject. 

The spiral spin state can be found in many physical 
systems. For example, both of neutron diffraction— and 
polarization dependent resonant soft x-ray magnetic scat- 
tering (RSXMS)— ■ experiments indicate an incommensu- 
rate superstructure with Q = (0.5, 0.1738, 0) at low tem- 
perature, where Q is the wave vector of spiral spins. In 
the direction of chain (b-axis) the spiral angle is about 
4>' = Q ■ b — 62°, which is the angle difference of two ad- 
joint spin. We will show our result in the case of <p' = 50° 
which is comparable to experimental finding. The dis- 
crepancy could come from inter-chain coupling. Though 
we also consider the inter-chain coupling in Sec|Vl its 
magnitude may not be practical. Further work is needed. 

Dzyaloshinskii-Moriya (DM) interactio n 23,24 not only 
plays an important role in the multiferroic material, but 
also produces many exotic physical phenomena. It is an- 
other mechanism which can give rise to spiral spin state. 
It can act as a vector potential on the spin wave in the 
magnon spin Hall effeclj— In ferromagnetic nanowires 
DM interaction has profound effect on the motion of do- 
main walls—. It can also give rise to spin current and 
soliton in spin chains— Therefore, it is important to in- 
corporate DM interaction into the model Hamiltonian to 
see what role it plays. 

The purpose of our study is to find the condition for 
the spiral spin state being the ground state (GS) in zigzag 
spin chain. In Scc|Hl we start with a Hamiltonian with 
NN, NN and DM interactions and derive the conditions 
of spiral spin states being the eigen states by performing 
a unitary transformation. In Sec lIIIl we use the pos- 
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itive semi-definite (PSD) matrix theorem to determine 
the conditions of ground state (GS), by decomposing the 
system into local Hamiltonian. Examples are given in 
Sec lIII Al and Sec lIII Bl to illustrate our result. The prob- 
lem of symmetry is discussed in Sec lIVI The applications 
to real physical systems of coupled zigzag spin chain and 
two-dimensional triangular lattice are given in SecJV] In 
Sec. IVI1 we compare our results with those of numerical 
calculations and simulations. Sec. I VIII is devoted to the 
conclusions. 
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II. SPIRAL SPIN STATE AS AN EIGEN STATE 

The physical systems which prefer spiral spin configu- 
ration usually have competing interactions. For example, 
when the nearest neighbor (NN) exchange interaction is 
weak, the next nearest neighbor (NNN) interaction or 
even Dzyaloshinskii-Moriya (DM) interaction become rel- 
atively significant. The Hamiltonian of this kind of sys- 
tems (are usually called the zigzag spin chains) can be 
written as 



H 



E 

3 



Ji[A lS | Sj 2 +1 -+ 
Dis-j x s j+ i ■ z] 



s j s j+i) 



J 2 [A 2 s?s| +2 + 
D 2 Sj x s j+2 ■ z 



^(sts j+2 + s jS j +2 ) 



(1) 



where j is the label of the lattice site, J\ is the NN in- 
teraction, J 2 is the NNN interaction, Ai(A 2 ) gives the 
anisotropic interaction along z axis for NN (NNN) inter- 
action, and Di(D 2 ) is the DM interaction between NN 
(NNN). 




FIG. 1: The way to analyze spiral state is to connect it to 
a fully polarized (FP) state with a unitary transformation, 
which is the product of rotations around z-axis. 

Although there have been numerous studies on the 
zigzag spin chain, the boundary for the spiral spin state 
being the ground state still cannot be determined if DM 
interaction is present. Here we propose another way to 
analyze. We connect the spiral spin state and a fully po- 
larized (FP) state by a unitary transformation. The FP 
state and the spiral spin state are shown in Fig. [TJ We 
then ask the question: Under what conditions the FP 
state will be the ground state (GS) of the transformed 
Hamiltonian? Since under any unitary transformation, 
the energy spectrum does not change, the ground states 
of the physical Hamiltonian and transformed Hamilto- 
nian are equivalent. By identifying the region of FP state 



FIG. 2: The definition of spin orientation n, where h 
(sin(/3) cos(a), sin(/3) sin(a), cos(/3)). 



being GS for the transformed Hamiltonian, we will know 
the region of spiral spin state as GS of the physical Hamil- 
tonian. 

The unitary transformation mentioned above rotates 
the spins around the z-axis^. It has the form 



N 



(2) 



with <f>' — Q ■ (Rj+i — Rj) being a constant. For now, 

(j)' and hence, Q can have arbitrary values. But, as it 
will be shown later, the DM interaction has a profound 
effect on <fi' and Q. For a given strength of DM interac- 
tion (Di and D 2 ), the boundaries of 0' and Q are deter- 
mined. Within the boundaries, the spiral spin states are 
the ground states of the system. Hence, the spiral spin 
state 

N 8 R - - 

\i>ss) = II[ cos (f )l t>j + S in(|)e i «^-*>| I),], 

3=1 

can be transformed into a FP state with spin direction 
n = (sin(/3) cos(a), sin(/3) sin(a), cos(/3)) (see Fig. [2]), as 
the following form 

N R R 

\FP,n) = [«*(§)! t)i +sin(|)e--| (3) 

The Hamiltonian is also transformed 
H rot = UHU- 1 



= £[J 1 (A 1 ^ +1 



cos(</>i — (j)') 
2 cos((/>i) 



( S ^ S J + l + S j S 3 + l) 



sm(0i -<?)_ _ 
H — - — Sj x Sj+i ■ z) 



cos(0i) 



cos(0 2 - 2</>') 



J2(A2S ^ + 2cos(0 2 ) W 8 7+* + S J S ^ ] 
sin(^ 2 - 20') 



cos ((fo) 



Sj X Sj + 2 ■ Z)] 
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where we have set <pi = tan 1 (Di), 2 — tan 1 (D 2 ), and 
used the identity 1 + itan(^) = sec(0)e J< ^. 

In order that \FP,n) is an eigen-state of the trans- 
formed Hamiltonian, it is required that 



ten as 



Ai = 
A 2 = 



cosf 



cos(0i) 
cos(<fc - 2<j>') 
cos(4> 2 ) 



(4) 
(5) 



The constraints give the limitation of the method, i.e., 
only under these conditions we can proceed further. On 
the other hand, the constraints show the relation between 
the spiral angle and the anisotropic exchange interaction. 
As a result of Eq. Q and Eq. J5]), H rot becomes an 
isotropic Hamiltonian 

H iso =2J JiAi[s*j • Sj + i + D^Sj x 8j+i ■ z] 

3 

+ J 2 A 2 [sj ■ Sj + 2 + D' 2 Sj x s j+2 ■ z] (6) 

where D[ = tan(0i - ft), D' 2 = tan(^ 2 - 2<f/). 

As the anisotropy exchange interaction is "rotated 
"away for both NN and NNN interaction, it is readily 
shown that \FP, n) is an eigen state of Hi SO with the 
relation 



x Sj ■ z\FP;h) = 0. 



To prove the above equation, one only has to rotate n to z 
and z to another direction. A more detailed analysis will 
be given in Sec. II VI where the symmetry of the system is 
also discussed. As a result, 



E a \FP,n) 



(7) 



where E = N(JiAi + J 2 A 2 )/4. In fact, Eqs. (4) and 
(5) combined is the requirement that NN and NNN ex- 
change interaction with anisotropy can be transformed 
into isotropic exchange interactions simultaneously. 



III. SPIRAL STATE AS THE GROUND STATE 



N 



Hiso - hjj+lj+2 



3 = 1 



JV-2 



= ^1,2,3 <8 l 2 N-3 + 22 123- 1 ® fljj+lj+2 <8> l 2 N-2-j 
3=2 

+ i 2 N-3 ® Hn-1,N,1 + i 2 f-3 ® /lJV,l,2 (8) 

where the local Hamiltonian is giving by 



hj,3+i 



3+2 



+ 



+ s j+2 ) ■ s j+ i] + J 2 A 2 Sj ■ s j+2 
JiAiD[ ^ _ „ 

[Sj X S j+1 + Sj+i x s j+2 \ ■ z 

J 2 A 2 D' 2 Sj x s j+2 ■ I (9) 



and 1m denotes the identity matrix of rank M. The di- 
rect product of hjj + ij +2 and 1m is meant to enlarge the 
vector space to 2 N dimension to accommodate TV spins. 
The direct sum gives a Hamiltonian matrix of 2 N dimen- 
sion, as shown in Fig. [3] If FP state of three spins is an 
eigen state of frjj+i j+2, then that of N spins is an eigen 
state of Hi SO , Furthermore, if a state of three spins has 
the lowest energy under Hjj+i j+2, then corresponding 
state (constructed by direct product) is the GS of Hi S0 . 
This is implied by the theorem of positive semi-definite 
matrix discussed below. 
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FIG. 3: The Hamiltonian can be decomposed into local 
Hamiltonians, and each local Hamiltonian contains only the 
interactions between three neighboring spins. 



In this section, we will identify the region of FP state 
as the GS in Hi SO (correspond to the spiral state as GS in 
H ) by decomposing the Hamiltonians into local Hamil- 
tonians and applying positive semi-definite theorem for 
analyzing. This method 29 had been applied to zigzag 
spin chains without DM interaction. It can also be ap- 
plied to spin of any length. Here we focus on spin-1/2 
system. 

We dissect the Hamiltonian into many local Hamiltoni- 
ans as shown in Fig. [3] Each local Hamiltonian contains 
three spins. Thus, the original Hamiltonian can be writ- 



To find the region for a FP state with spin direction 
n as the ground state of Eq. (|5J|, we rotate the 2-axis 
to direction h' = (sin(/3) cos(— a), sin(/3) sin(— a), cos(/3)) 
and n to z. The Hamiltonian of Hi SO in Eq. (|9]) becomes 

h j'j+l,j+2 = ~ ^K^' + ? i+ 2 ) ' Sj+l] + J2A2S3 ■ S j+2 

JiAiD[ „ _ ~, 

+ ^ ^' X + S 3 + 1 X S l+ 2 \ ' H 

+ J 2 A 2 D' 2 Sj x Sj + 2 ■ n'. 
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and 



^ J + l J +2lt> J lt> J + l|t>,+2 

= S 2 (J!A 1+ J 2 A 2 )|t>,|t>, + i|t) 



J+2 



^sin (/3 )e-(^i + J 2 A 2 ^) 



4i " J 2 

x (It^ltWU^-U^lt^+ilt)^) 



where = tan(0i — 
found that for | f)j\ 
the required relation is 

J1A1 
J 2 A 2 



tan(0 2 - 20'). It is 



i| t)i+2 to be an eigen state, 



-2J> 2 



(10) 



This procedure is not restricted to spin one-half systems. 
Substituting Eq. (flU)) into Eq. ([5]) , the local Hamiltonian 

becomes 



D 1 

'ljJ+l.j+2 = J2A 2 [-^ r (Sj + S j+2 ) ■ S j+ i + Sj ■ S j+2 



- D' 2 (sj x Sj + i + Sj+i x s j+2 ) ■ z 
+ D' 2 Sj X s j+2 ■ z] (11) 

To find the region for FP state to be the ground state 
of this local Hamiltonian, we applied the positive semi- 
definite theorem. The theorem of positive semi-definite 
(PSD) matrix is described by the following: The neces- 
sary and sufficient condition for a real symmetric matrix 
A to be positive semi-definite is x T Ax > for all real 
vectors x. If M and N are positive semi- definite, then 
the sum M + N, the direct sum M© JV, and direct product 
M ® N are also positive semi- definite. Hence, if we are 
able to prove that (/ij j+i,i+2 — Eq) is a PSD matrix with 
Eq being the energy of the FP state, then the Hi SO — Eq 
is also a PSD matrix. For spin 1/2 system, the Hamilto- 
nian in Eq. (fTTj) is 2 3 x 2 3 matrix. Its energy spectrum 
is 

Eq with four fold degeneracy; 
Eq + 5 Ei with two fold degeneracy 
Eq + SE 2 with two fold degeneracy (12) 

where Eq is the energy of FP state, and 



6E 1 _ 2D[D' 2 - D' 2 - ^[2,D' 2 D 12 + (D[ + D 12 )] D' 2 



J 2 A 2 



2D' 2 



SE 2 _ 2D[D> 2 - D' 2 + VKfl? + {D' 1 +D' 2 )} D' 2 



J 2 A 2 



2D' 2 



To make local Hamiltonian a PSD matrix, one requires 
6 Ei > and 5E 2 > 0. Therefore the conditions for PSD 



This is the end of our derivation. Summarizing briefly, 
Eq. (4) and Eq. (5) are the conditions of FP states being 
the eigen states of Hi OS in Eq. (3). Hence they are also 
the conditions of the spiral spin states being the eigen 
states of the physical Hamiltonian in Eq. (1). On the 
other hand, Eq. (10) and inequality (13) are the condi- 
tions of the spiral spin states being the ground states of 
the physical Hamiltonian in Eq. (1). 

For most insulating compounds, J 2 , the NNN superex- 
change interaction is antiferromagnetic. Hence, we con- 
sider the case J 2 > and A 2 > 0. The region for PSD 
in Eq. (fl~3|) can also be written as 



2(f)' - f < 2 < 20' 
2<t>' < 02 < 20' + f 



for S2. < 



for 



< 



< 
< 



<l>2 

2 



It is indicated by the shaded area in Fig. [4j for a given 
0'. Fig. |4] shows the main result of this work. In view of 
Eq. flU, Eq. ([5]) and Eq. (flQ|) . we see that there are two 
free parameters in Eq. (TTJ), namely 0i(Z?i) and 02(-C>2)- 
In the parameter space of 0i and 2 , the shaded regions 
show where spiral spin states are the ground states. The 
pitch angle 0' of the spiral spins is closely related to the 
values of 0i(Z?i) and 2 (I?2)- Hence, the DM interac- 
tion plays the crucial role of determining the existence of 
spiral spin states. 

It has been shown that inside the shaded region of 
Fig. the spiral spin states are the ground state. How- 
ever, we also found that there are gapless excitations or 
Goldstone modes (details are in sec. ITV)) . Outside the 
boundary, these modes actually have lower energy as it 
will show in Eq. fp~9|) . Hence, the spiral spin states are 
the ground states in the shaded region. But they are not 
the ground states outside of these boundaries which are 
exact since our solution is exact. On the other hand, it 
is possible that the chiral correlation or the in-plane spin 
correlation still exhibit long-range order behavior out- 
side of the boundaries. Furthermore, we cannot rule out 
the possibility of the spiral spin states being the ground 
states in some other region in phase space which is not 
in the neighborhood of the shaded regions. 

For later use, we will show some cases explicitly in Fig. 
O the shaded regions will correspond to the condition for 
spiral GS. In the following, we give two simple examples 
to illustrate our result. 

Similarly, for spin-s system of the Hamiltonian, the lo- 
cal site Hamiltonian is a (2s + l) 3 x (2s + l) 3 matrix. 
After diagonalization, the energy spectrum can be ob- 
tained. The theorem of PSD matrix can also give the 
region for FP state as GS. Hence the result of this sec- 
tion can be generalized to spins of any possible length. 



(1 - D' 2 )D' 2 > 2D[D' 2 and 2D[D' 2 > D' 2 for J 2 A 2 > 



A. Isotropic exchange interaction 



(1 - D' 2 )D 2 < 2D[D' 2 and 2D[D' 2 < D' 2 for J 2 A 2 < 0. For physical systems having isotropic superexchange, 

(13) or Ai = A 2 = 1 (in our case 0i = \,4>2 — <t>'), the 
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0'+f 



<t>' 



<t>'- 




2<p' - 



2o' 



FIG. 4: The shaded regions show where spiral spin states are 
the ground states and the relation between the pitch angle 
and (0i, 02) in the Hamiltonian, and the region for shifting n 
in both (0i and 02) axes are also the PSD region. 



Hamiltonian of Eq.([T]) becomes 

, v-^ tan(0') 
n = <h ~ 2 : jjr^ s i ■ s i+1 + Sj • s i+2 



tan(- 2 , 

tan(0')(2s*i x s i+ i ~~ Si x s l+2 ) ■ z 



The conditions for PSD of Eq . (]13|) turn out to be 
0° - 90°, 270° - 360° with J 2 > 0, 



90 u - 270 u 



with J 2 < 0. 



(14) 



(a) 0' = 3O° 



2tt 
3 




FIG. 5: The shaded regions show where spiral spin states are 
the ground states and the relation between the pitch angle 
and (0i, 02) in the Hamiltonian, for (a)0' = 30°, (b)0' = 50°. 



B. Without DM interaction 



So in physical case J 2 > 0, the spiral angle 0' is in the 
range 0° - 90° or 270° ~ 360°. 

The implication of Eq. (fT4"]) can be seen by con- 
sidering the following three simple cases: (a) 0' = 
0°, or the GS is the FP state. In this case H' = 
J 2 J2i ~4si • s*j_|_i + Si ■ Si + 2- It has been showr*i£ that 
the FP state and the UDRVB state are degenerate 
GS. Our method correctly leads to one of them, (b) 
0' = 7r, or the GS is the Neel state. In this case 
H' = — \J2\^iSi ■ Si+2- The system becomes two de- 
coupled spin chains. Each is ferromagnetically coupled, 
(c) The GS is spiral state with spiral angle being 30°. 
The resulting parameters are ( Ji, J 2 , A1, A 2 , D x , D 2 ) = 
(-4.13,1,1,10.268,0.577). (d) For specific spiral 
angle 0' = 50°, the resulting parameters are 
(Ji, J 2 ,Ai,A 2 ,£>i,£> 2 ) = (-5.111,1,1,1,0.466,1.191). 

In this isotropic Hamiltonian, the relative strength of 
exchange coupling is J1/J2 = — 2tan(0')/tan(0'/2). In 
the region of PSD in Eq. (TT4")) with J 2 > 0, this is always 
less than —4. We shall see in Sec. [V] that for coupled 
spin chains or two-dimensional cases, the superexchange 
Ji I J 2 can be tuned by inter chain coupling J3 to make 
|Ji/J 2 | smaller. 



For system containing high symmetry or small spin- 
orbital interaction, we have the case D\ « D 2 w 0, im- 
plying 0i = 2 = 0. Then the Hamiltonian in Eq. ([T]) 
under the condition of Eqs. ()4I5|10[) becomes 

H " = j 2 Y^ -4cos(0')(cos(0') s ^ +1 + s * s * +1 + s y s y +1 ) 

3 

+ (cos(20>!4 +2 + s* s * +2 + s y s y +2 ) (15) 

and the spiral state is the eigen state. Now we have 

^=-4cos(0'), (16) 

which is consist with the result of treating spins as classic 
vectors. 

Here are several examples of the system 
in Eq. (fT5l) : (a) For GS with pitch angle 
0' = 30°, the Hamiltonian has parameters 
(J 1 ,J2,A 1 ,A 2 ,D 1 ,D 2 ) = (-3.464,1,0.866,0.5,0,0). 
(b) For GS with specific pitch angle 0' = 50°, the 
Hamiltonian has parameter (Ji, J 2 , Ai, A 2 , D%, D 2 ) = 
(-2.571,1,0.643,-0.174,0,0). 

The spiral spin state is always the ground state pro- 
vided the anisotropy interactions between NN and NNN 
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have the forms in Eq. ([15]) with J 2 > 0. However, 
this kind of solution may not be possible for multi- 
ferroics LiCu2C>2 with <\>' = 62°, because in this case 
A-2 = cos(124°) < 0, and it is likely to be unpractical. It 
will be shown later, even if inter-chain coupling is con- 
sidered, A2 remains negative. This suggest that there is 
DM interaction in LiCu202- 



C. Special case with 4>' = 50° 

For multiferroic compound LiCri202, the spiral angle is 
close to 7r/3, therefore here we consider this special case. 
DM interaction and anisotropy are likely to be small. 
The parameter set {(J x , J 2 ), (A x , A 2 ), (D x , D 2 )} = 
{(-3.111,1), (1,0), (0.176, 0.176)} gives us 0' = 50°, 
which is close to ir/3. However, we note that A 2 vanishes 
in this case. We suggest that one can treat J2^2S z s z +1 
as a perturbation. Since J2A 2 > 0, the perturbed state 
should have greater spiral angle. Hence, our approach 
should to be capable of dealing with realistic physical 
system such as LiCu202. The same approach can be 
used for other value of spiral angle </>'. 



IV. SYMMETRY PROPERTIES 

In Sec. mi we showed that FP state in Eq. ([7]) for any 
direction n is an eigen state. Later in Sec. IIII1 we iden- 
tified the region for these states to be GS. These states 
form a Hilbert space which has SU(2) symmetry. We 
can implement a site-dependent unitary rotation (0(2) 
rotation) to generate spiral states. On the other hand, 
the Hamiltonians H, H rot and Hi SO all have only SO(2) 
symmetry. Hence, we arrived at a situation which is 
called "emergent symmetry" by Batista 30 -. The symme- 
try group of the Hamiltonian in Eq. ([1} is isomorphic to 
that of Hi SO which is 50(2). However, the Hilbert space 
of degenerate spiral states has a symmetry group which 
is isomorphic to SU(2). We give a detailed analysis in 
this section. 

To see the symmetry property of this system more 
clearly, we first express our results following the nota- 
tions of Batista 3 - -. We note in passing that our result 
contain DM interaction, which is a generalization. The 
Hamiltonian of Eq. ([1} in fc-space is 

H = ^[JiAi cos(q) + J 2 A 2 cos(2q)]s z q s z _ q 
<? 

+ [Ji(cos(g) - D 1 sin(g)) 

+ J 2 (cos(2g) - D 2 sin(2<?))]s+s- (17) 

where S * = ^ e«^, •+ = ^ ^sj, and s q = 
(s'g'y , with lattice constant set to 1. The commutation 
relation of s q=< p, is 



where 



t,H] 



7z 



ai 



Jisin(0i - (/>') 



(s 



1+1 



J 2 sin(<£ 2 - 2(f>') 



cos 



{8 



1+2 



-l) 



"1-21 



From H\FP, -z) = E \FP, -z) and m\FP, - 
get the relation of H in Eq. (fT7l) as 

[s+,(H-E )](s+y\FP,-z)=Q 



0, we 



for any integer p. Hence, (sJ,) p |FP, — z) is one of the 
degenerate eigen states of H . The spiral spin states is 
a linear combination of (sj, ) P \FP, z), we show it more 
directly below. 

It is clearer to discuss the symmetry property by ana- 
lyze the transformed Hamiltonian. The FP states are the 
degenerate GS of Hi SO irrespective of the direction of n. 
To put it in another way, H iso commute with J2i s t so 
that the z-component of the total spin is a good quantum 
number, hence, they have the following eigen states 



\ 3 ,k,l),...,\FP;z)} (18) 



j<k 



j<k<l 



where \j,k,...) = S+S£...\0) and |0) = \FP;-z). 
Since any FP state \FP, h) is a linear combination of 
the above states, the states in Eq. (fT8|) are actually de- 
generate GS of Hi SO , Under the unitary transformation, 
the above states can be transformed into magnon states: 

{|o),4|o),(4) 2 |o),(4) 3 |o),...,|pp;^)} 



The states (s^,) p \FP,-z}, those in Eq.$TE§ and \FP;h) 
are three basis of the degenerate GS of the Hamiltonian in 
Eq. (JTJ) . Each set form a Hilbert space with a symmetry 
group isomorphic to SU(2) which contains the symmetry 
group (isomorphic to SO(2)) of the Hamiltonian. 

As for the excitation energy of the Hamiltonian in Eq. 
([T]), consider one magnon with wave vector k as \ip) = 
Tlj e Using the relation 

(H - NE )\j) = -(JxAi + J 2 A 2 )|j) 

+ ^.[( 1 + iD l )\j-l) + (l-iD 1 )\j + l)} 

+ ^[(l + i J D 2 )|i-2) + (l-i£» 2 )|i + 2)] 

we found that |^) is an eigen state and its energy spec- 
trum is 

w k = E - NE 

= Ji(cos(fc) - Ai + Di sin(fc)) 

+ J 2 (cos(2/c) - A 2 + L> 2 sin(2fc)). (19) 

When k approaches </>' , the energy E converges to Eq. 
Hence we have gapless excitation. Note that the one- 
magnon states are not the only low-lying excitations. 
In fact, there are certain two- magnon, multiple- magnon 
states that are also gapless excitation. Studying on this 
topic will be presented in another paper. 
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V. EXTENSION OF THE MODEL 

Our result can be generalized to the cases of coupled 
spin chains and certain types of two-dimensional trian- 
gular lattices. Hence, physical systems such as LiCu2C>2 
and ACr02(A=Cu,Ag,Li or Na) can be studied with our 
method. 



A. Coupled spin-1/2 zigzag spin system 



Ji Hi J3 




(a) GO 



Ai(A2,A3) is the anisotropic interaction along z axis 
for NN (NNN, inter-chain) interaction and D^s are the 
respective strength of DM interaction. Superscript esc 
denotes coupled spin chain. With similar procedure 
in SecHH we set D\ — tan(0i),D2 = tan(02),Z?3 = 
tan(03), and construct the unitary transformation rotat- 
ing the spins around z-axis: U — Ylf—i exp(is|Q • Rj) 

with a constant 0' = 2Q ■ (Rj+i — Rj). The spiral angle 
is shown in Fig. [SJb). This gives a fixed phase difference 
between two chains. The conditions of finding a uni- 
tary transformation /7(0') to change H into an isotropic 
Hamiltonian are 



cos(0i 


-00 


cos( 


h) ' 


COs(0 2 


-20') 


cos( 




COS(03 


2 > 


cos( 





The resulting isotropic Hamiltonian is 

i 

+ J 2 A 2 [Sj ■ Sj+4 + D' 2 Sj X Sj+4 • z] 

+ J3A 3 [sj ■ Sj+i + D' 3 Sj x s j+ i ■ z] 



FIG. 6: (a) Two coupled zigzag spin chains, (b) Two basic 
elements denoted as a blue (light) triangular and a red (dark) 
triangular for two coupled spin chains. 

Recently, the compound LiCu 2 2 11112 attracted re- 
searchers' attention because of its multiferroic property. 
It has chains of edge-sharing oxygen plaquettes with cop- 
per ions at the centers. The NN and NNN superexchange 
interactions strength between copper ions are compara- 
ble. A coupled zigzag spin ladder with two legs are shown 
in Fig. [nja). The model Hamiltonian is given by 

HCSC =EM*is*s* +2 + \{sp- +2 + sJsf +2 ) 
j 

+ D 1 Sj X Sj +2 ■ z] 

+ J 2 [A 2S | S | +4 + ^(4*7+4 + 574+4) 
+ D 2 Sj x s j+4: ■ z] 

+ D^Sj x Sj + i ■ z] (20) 

where j is the label of the lattice site, J\ is the nearest 
neighbor (NN) interaction, J 2 the next nearest neigh- 
bor (NNN) interaction, J3 is the inter-chain coupling, 



where D[ = tan(0i - 0'), D' 2 = tan(0 2 - 20'), and D' 3 = 

tan(03 — \)- If for Hamiltonian H£££, FP states are the 
eigen states then for Hamiltonian H csc , spiral states are 
the eigen states. 




h,A h B 



FIG. 7: The two basic elements of the system:red (dark) tri- 
angles Ka and blue (light) triangles hs- 



We again dissect Hf^ into local Hamiltonians, each 
contains three spins. However, now there are two kinds 
of local Hamiltonians, represented by blue (light) and 
red (dark) triangles in Fig. EJb) . Since the two triangles 
have a common edge, we divide J\ into J\ a and Ju- This 
division is arbitrary. As shown in Fig. the two local 
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Hamiltonians are Ha (dark or red)and hs (light or blue:) 

flA = ^ o 1 [(Si + S i+i ) ■ S i+2 ] + J2A 2 Si ■ s l+i 



JiaAiDi 



J 2 A 2 D 2 Si X S i+4 • £ 

J3A3 



(21) 



2 

J3A3D3 



[(s*i + S i+2 ) • S i+1 ] + JlfeAxS, • S i+2 

s« x s i+ i + si+i x s i+2 ] • 5 



JlbAiD^S,; X S i+2 • 2. 



(22) 



l-FP, n) will be the eigen state locally, provided the follow- 
ing conditions which are analogous to that in Eq. (ITU1) . 
are satisfied by 



(23) 
(24) 



Jla 


-2D' 2 A 2 


jT ~ 


D[Ax 


A- 


-2£>iAi 


Jib 


D' 3 A 3 



The forms of both local Hamiltonians Ha and hs in Eq. 
([2T]l and Eq. ((22j) are the same as those of Eq. ([IT]). 
Hence the spectrum in Eq. (fi"2"]) is applicable. The energy 
difference between the FP state and others for Ka and hs 
are 



<5£± _ 2£>i£> 2 - Df ± ^[3D^ 2 L>i 2 + (L>i + D0 2 ]Df 



J 2 A 2 
JibAi 



2L>i- 



/2 



<5£± 2£>^ - D' 2 ± ^/[3D' 2 D' 2 + (D' 3 + D'^Df 

2Df 



For H — Eq to be a PSD matrix, one requires that 5E^ > 
and 5E„ > 0. Therefore the conditions for PSD are 



and 



(1 - D' 2 )D' 2 > 2D[D' 2 and 2D' 2 > D[ for J 2 A 2 > 0, 

(1 - D' 2 )D' 2 < 2D[D 2 and 2D' 2 < D[ for J 2 A 2 < 

(1 - Df)D[ 2 > 2D' 3 D[ and 2D[ > D' 3 for J 16 Ai > 0, 

(1 - D I2 )D 12 < 2D' 3 D[ and 2D[ < D' 3 for J lb A x < 0. 



The implication of above derivation can be seen by the 
cases similar to those of Sec lIII Al Scc lIII Bl and Sec lIII Cl 

• For the isotropic case with Ai = A 2 = A3 = 1, and 
hence 4>i = 4>'/2,(f> 2 = 4>',(f> 3 = 0'/4, the Hamilto- 
nian is 

HCSC = Et— 7^( J 2 tan(0') + J 3 tan A]*! ■ s t+2 
i tan( — ) 

+ J-lSi ■ Si+4 + J 3 Si ■ S i+ i 

+ 2[J 2 tan(0') + J 3 tan(— )]sj x s l+2 ■ z 
+ J 2 tan(<//)s*i x s i+4 ■ z 

+ J 3 tan( — )Si x s i+ i ■ z 

The resulting relations between J\ and J 2 , J3 from 
Eq. j23j and Eq. ([24]) are 



Ji 



Similar to Eq. ([T4"]) . we get the PSD conditions for 
h A in Eq. j2T} and h B in Eq. ([22]) 



and 



0° - 90°, 270° - 360° with J 2 > 
90° - 270° with J 2 < 

0° - 180°, 540° - 720° with J lb > 
180° ~ 540° with Ji b < 



— — [J 2 tan^') + J3tan(^-)]. (25) 



Note that Q • (R i+ i - R t ) = (f>'/2, so 0' and 
0' + 27r do not give the same spiral spin state. 
For antiferromagnetic coupling (J 2 > 0), and 
no constraint on Jib, the PSD condition re- 
quires the spiral angle <f>' to be in the region 
0° - 90° or 270° - 360°. We give the follow- 
ing examples. For spiral angle (jJ — 30°, and 
J3/ Jn = ±0.3, the following set of parameters sat- 
isfy PSD condition: {(J 1 ,J 2 ,J 3 ),(D 1 ,D 2 ,D 3 )} = 
{(-3.772, 1,-0.3), (0.268, 0.577, 0.132)} with posi- 
tive J 2 and Jib- 

For the case of no DM interaction as that in 
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Sec. IIIIB1 we again give the example of 4>' = 
30°, J3/J2 = ±0.3. The parameters satisfying PSD 
conditions can be {(Ji, J 2 , J3), (Ai, A 2 , A 3 )} = 
{(3.386, 1,-0.3), (0.866, 0.5, 0.966)} with positive 
J 2 and Jib- 

For the case $ — 50°, we set <j>\ = fa = fa = 10° 
and J3/J2 = ±0.3, and get the parameter set of the 
system {(J 1} J 2 , J 3 ), (A 1; A 2 , A 3 ), (D u D 2 , D 3 )} as 
{(-3.151, 1, -0.3), (0.778, 0, 1), (0.176, 0.176, 0.176)} 
with positive J 2 and Jib- Above parameters are 
very close to the realistic systems, For example, in 
LiCu 2 02, the exchange interactions given by ab 
initio calculation^! are very close to our parame- 
ters. The DM interaction we get are smaller than 
J 2 by an order of magnitude. The only possible 
discrepancy is A 2 which should be close to 1. We 
suggest the term J 2 A 2 sf sf +2 should be treated as 
a perturbation. 



B. 2D triangular spin-1/2 spin system with spiral 
state 




case J3 = J 1 , then we have 

H 2D 



j,k S j',k 



s j,k s j',k) 



+ J 2 [A 2 s| )it sf >fe , + -(s+ 
+ D 2 s^ k x s,j^ ■ z\] 



],k S j,k' 



S j,k S t,k' 



(26) 



where j stands for the site index along 6-direction, 
and k is the index along a-direction. Defining Di = 
tan(^i), D 2 = ta,n(fa), and performing a unitary trans- 
formation U(<f>'), we get an isotropic Hamiltonian 



H 



2D 



(j,k) 



with the parameters Ai = 
and D[ = tan(fa - </>'), D' 2 



J 2 A 2 [sj,fe • Sj,k> + D' 2 s-j : k x Sj,k> ■ z] 

cos(4> 1 —4>') ^ cos(0 2 — 20') 







COS(0i) ' 

tan(^ 2 - 24>'). For Hamil- 
tonian HfJZ , FP state is the eigen state, and for Hamil- 
tonian H 2D , spiral state is the eigen state. 

The Hamiltonian of a layer is the combination of many 
local Hamiltonians of small triangles, as shown in Fig. 
[8jb). For each local Hamiltonian, we can perform our 
analysis as before. The local site Hamiltonian is the same 
as that described in ID zigzag spin chain. Therefore, the 
region for PSD is the same as that of zigzag spin chain 
in Eq. (fT3"|) . and Eq. (fit))) also applies. 

Consider isotropic exchange interaction as an example. 
We set A x = A 2 = 1 and $1 = 0'/2,</> 2 = <t>' ■ It was 
found that the PSD condition gives 



U 



90°, 270° ~ 360° 
90° - 270° 



for J 2 > 
for J 2 < 



, .2. 



Our method has the potential application to two- 
dimensional spin system. 



7/ \ 2 



2 



Js. \ /A 

2 \ '2 



(b) 



FIG. 8: (a) Two-dimensional triangular lattice system con- 
taining Ji(red), J3 (black) and J2(blue) superexchange inter- 
action, with the spiral angle at different sites are shown, (b) 
Two basic elements of the triangular lattice. 



Our method can be applied to the compounds of layers 
of triangular lattice, such as ACr02^(A=Cu, Ag, Li or 
Na) shown in Fig. Ufa). Here wc assume that the ex- 
change interaction Ji along a-axis,J 2 along 6-axis is J 2 
and J3 along a — b axis. For simplicity, we consider the 



VI. ANALYSIS AND DISCUSSION 

We now compare our results with those of numeri- 
cal calculations and simulations. The results from the 
calculations (see below) show that the zigzag spin chain 
possesses rather rich and sophisticate physics. Since our 
method is to find the exact solutions, we can only access 
limited regions in phase space. However, inside these re- 
gions, our results are exact and thus can be used as stan- 
dard to check the results of numerical calculations. Our 
calculation can also serve as a guideline for the calcula- 
tion of the correlation functions outside of these regions. 
Hence, by comparison, we hope to shed some light on the 
complex physical landscape of the zigzag spin chain. 

There are many works on the phase diagram. Hiki- 
hara et ali^ used Density matrix renormalization group 
(DMRG) method and found spin-liquid, dimer and gap- 
less chiral phase for various quantum spins. For s — 
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1/2, they found notably, the chiral correlation function 
(sj x sj +( i), has long-range order behavior. Plekhanov 
et al.— also used DMRG technique to study similar sys- 
tems. They found that quantum fluctuation modified the 
classical phase diagram. In particular, they found phases 
denoted by spin liquid I, spin liquid II, E-I and E-II. In 
the former two phases, they found that the spin corre- 
lation functions (especially those on x-y plane (sjsj +d )) 
have the power-law decay behavior. In the E-I and E-II 
phases, the spin correlation functions decrease exponen- 
tially except for a weak magnetization in E-I phase. Ja- 
fari et ali^i used quantum renormalization group, to find 
a similar phase diagram except for a new phase near the 
spin liquid I. It was designated as dimer II. For J\ < 0, 
the spin liquid III have more or less occupied regions of 
E-II phase of Plekhanov et alj^ 3 .. Dmitriev and Krivncv*^ 
studied the same problem with variational mean-field 
approximation and scaling. They also obtained simi- 
lar region of spin-liquid I in the phase diagram. On the 
other hand, their incommensurate phase (III) fits roughly 
within the region of E-II phase of Plekhanov. 

Since all the above calculation did not consider DM 
interaction, we set D\ = Z? 2 = at beginning in order 
to make comparison. This, in our calculation, is to set 
Ai = cos(0'),A 2 = cos(2</>'). The exact solution can 
be found when Eq. (fT6|) is satisfied. Wc consider the 
case J2 > 0. The notations of the phase diagram of 
Dmitriev and Krivnov^ was chosen because their mean 
field calculation gives clear boundaries between phases. 
Besides, the result of Plekhanov et al£& is consistent with 
that of Dmitriev and Krivnov^. Our result is shown by 
the heavy solid line in Fig. |H1 The ground states are 
the spiral spin states with different pitch angles. The 
exact solution is in the region of incommensurate phase 
or the E-II phase of ref. 33 or that of spin liquid III of 
ref. 35. Our calculation shows that on the solid line the 
one-magnon state is gapless and the chiral correlation has 
long-range order behavior due to spiral spins. For larger 
Ai, numerical calculations, except for ref. 33, found a 
gapless chiral phase. For smaller Ai , it has been argued 34 
that it is a phase (denoted by spin fluid I) of XXZ model 
without long-range order. Our calculation also shows 
that one-magnon modes actually have lower energy than 
that of spiral spin state as it can be seen from Eq. (TT91 . 
Hence, it is plausible to assume that the solid curve is the 
exact boundary of a gapless chiral phase. We noticed that 
our curve is close to those given in ref. 35 in the region 
0.25 < J = — J 2 /Ji < 0.28 but the difference increases 
as J2/J1 increases. Since the calculation of Dmitriev and 
Krivnov 3 -^ is less accurate for larger J 2 / Ji whereas ours 
is exact, this discrepancy is expected. Further study is 
needed to confirm that the curve by our calculation is 
indeed the phase boundary. A final note of this case is 
that there exists another "mirror" curve if one makes the 
transformation: Ai — > — Ai, J\ — > —J\,fa — > ir — fa. 
This is actually equivalent to rotating the spins of odd 
sites around z-axis by n. 

Next, we discuss the effect of DM interaction. With 



DM interaction, the general form of the ratio of the cou- 
plings changes from Eq. (flTj]) to the following 

Ji = 2 sin(<fe - 2fa) cosjfa) 
J 2 &m((f>i - fa) cos(</> 2 )' 

as was given in Eq. (|10j) . The system is more apt to have 
a spiral spin ground state. For example, in Fig. |H1 the 
heavy dashed (Di = 0.1) and dotted lines (D2 = 0.1) 
indicate the conditions of such ground states. In these 
calculations there is a small but non- vanishing A 2 . For- 
tunately the result of Dmitriev and Krivnov 35 shows that 
the effect of A 2 will not influence the phase diagram very 
much. The curves evidence that the spiral spin states can 
be the ground state in the spin-fluid I phase of ref. 33. 
In other words, DM interaction makes chiral correlation 
more robust. 

For the case J\ > which was considered by Hikihara 
et al<22, if there is no DM interaction, there is no spiral 
spin ground state. The reason is simple. In view of Eq. 
(fTBl cos(0') has to be negative for a positive J\j 'J 2 . This 
implies that Ai = cos(0') < 0, is contradict to our start- 
ing point of Ji > 0. But if the DM interaction is present, 
the spiral spin states can still be the ground state as long 
as the condition Eq. (|27|) is satisfied. For example, we 
can choose fa . The region of PSD is shown in FigJSJ We 
may have many possibilities. Some of them are 

(1.244, 3.50, 0,3.73, 0.176), 0i = 75° and 4> 2 = 10°. 
(5.85, 1.97, 0, 1.73, 0.176), fa = 60° and fa = 10°. 
(1.28, 3.50, 0.185, 3.73, 0.364), fa = 75° and fa = 20° 

where the parentheses give the values of 
{J1/J2, Ai, A 2 , Di, D 2 ). Under these conditions, 
we have J\ > 0, J 2 > 0. Hence, due to the effect of DM 
interaction, the spiral spin states can be the ground state 
in the regions of the spin-fluid phase and dimer phase 
designated by Hikihara et al.— even if the magnitude of 
the DM interaction should not be small. 

There have been many numerical calculations on the 
zigzag spin chain system. However, due to the large num- 
ber of parameters, and hence, large dimension of phase 
space, and the different methods of calculation, it is not 
always clear how physical properties change due to the in- 
terplays among these parameters. For example, how the 
system evolves if Ai A 2 < but its magnitude increases? 
A very helpful way is to calculate the correlation func- 
tions as it has been done in many works. In many cases, 
finding exact solutions, albeit in a very limited region 
can provide more clues. A simple calculation shows that 
the spiral spin state is able to sustain not only the vector 
chiral spin order but also the oscillatory spin correlation 
on the easy-plane 

<X>K+rf> = Y sin2(/3)cos(<//) ' 

j 

where fa is the spiral angle and j3 is the cone angle. The 
correlation should oscillate but not decay. Numerical 
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analysis should be able to find long-range order behavior 
of the chiral correlation function and the xy-plane oscil- 
latory spin correlation function. But there is some dif- 
ficulty due to the symmetry of the Hilbert space caused 
by the degenerate ground states as we have discussed in 
section IV. As it has been mentioned, if the anisotropic 
exchange interaction can be transformed away then it can 
be shown that all the FP states are the ground states 
of the transformed Hamiltonian. This implies that our 
Hamiltonian has huge degeneracy. The ground state can 
be any linear combination of these degenerate states. If 
in numerical calculation, the sample of ground states is 
not large enough, the symmetry can be lost and correla- 
tion functions decay faster than they should be. 

A=Aj 



Ref 33: A 2 =1 (No DM) 




Ref 33: A-|=A 2 (No DM) 




Ref 33: Classical result (No DM) 


Ref 33: MFA (No DM) 




Dl=0.1, D 2 = 






• ••• D-|=0, D 2 =0.1 




FM * 


No DM 




1 C 




SF (1) \ \ — - 

\ V N 

\ V 





0.0 0.2 0.4 0.6 0.8 

J = - J 2 / J 1 



VII. CONCLUSIONS 

To find the conditions for spiral spin state as the 
ground state, we transform the physical Hamiltonian to 
a Hamiltonian with isotropic exchange interaction and 
DM interaction. The FP states are the eigen states of the 
transformed Hamiltonian. Then we used positive semi- 
definite theorem to identify the region of FP state being 
the GS for the transformed Hamiltonian, which is noth- 
ing but the same region of spiral spin state as GS of the 
original Hamiltonian. The region (shown in Fig. |4j can 
be expressed by a very simple relation with the couplings 
of NN and NNN superexchange interaction and DM in- 
teraction. The effect of DM interaction is important be- 
cause its strength is related to the pitch angle of spiral 
spins and the unitary transformation. As the strength of 
the DM interaction increases, the pitch angle 0' also in- 
creases. It was also found that for spiral spin states to be 
GS, either 20! = 2tan~ 1 (D 1 ) > 20 2 = 2tan" 1 (D 2 ) and 
01 < 0' < 02/2 + 7r/4, or 2 /2 - tt/4 < 0' < 2 /2 and 
0' < 0i < 02/2 + 7r/2. Wherever the equal signs stand, 
the spiral spin states are degenerate with one-magnon 
states for the physical Hamiltonian in Eq. ([1]). Hence, 
the boundary in Fig. 3] marks the region where the spiral 
spin states being the GS. We thus made manifestly the 
relation between DM interaction and spiral spin states. 

These results can serve as a guided line for experimen- 
talists to find spiral spin state which, in turn, can lead to 
multiferroics. Our results also show the connection be- 
tween spiral spins and magnetic frustration. Finally, our 
method can be applied to other types of magnetic system 
such as coupled spin chains and layer structure. It is not 
restricted to spin- 1/2 system, but can be applied to any 
spin system. 



FIG. 9: The phase diagram according to ref. 35. Our results 
are shown in the solid line (no DM interaction), the dashed 
line (Di = 0.1, D 2 = 0), and the dotted line (Di = 0, D 2 = 
0.1). SF (I) and IC denote spin-fluid (I) and incommensurate 
phases. 
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